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Abstract 


In this report we study two problems related to the effects of 
turbulence in the atmosphere of Venus on the Pioneer Venus entry probe 
radio link. In the first, the cross correlation between the log~amplitude 
and phase fluctuations of the Pioneer Venus communications link is examined. 

It is seen that for fluctuation frequencies above approximately 1 Hz there 
is little or no correlation. For frequencies below, the correlation is 
weak and the square root of the coherence has a peak value close to 0.65. 

The second problem consists of inferring the turbulence characteristics 
of the Venus atmosphere from the Mariner 5 phase fluctuations. It is seen 
that with the data processing techniques developed and currently available, 
the phase error due to oscillator drift, assumed trajectory delay and 
spline-curve fit exceed the turbulence induced fluctuations. It is, therefore, 
not possible to infer the turbulence characteristics from the Mariner 5 
phase fluctuations. 
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I, Cross Correlation Between the Amplitude and Phase Fluctuations 


of the Pioneer Venus Entry Probe Communications Link 


1.1 Introduction 

Before the effects of a turbulent medium on a communications link can be 
realistically assessed, the statistics and characteristics of both amplitude and 
phase fluctuations induced in a radio signal propagating through the turbulent 
medium must be known. In our previous study (Ref. 1), we discussed in detail 
the distributions, means, variances and power spectra of the amplitude and 
phase fluctuations in the Pioneer Venus entry probe radio link. Also of 
interest is the degree to which the amplitude and phase fluctuations are 
correlated (Ref. 2); in other words, we wish to obtain information on the 
cross spectrum and coherence between the amplitude and phase fluctuations. 

Although the formulation for analyzing these correlation characteristics 
has recently been given (Ref. 3), no detailed results have been computed or 
published. The purpose of this report is to present calculations carried 
out and applied to the Pioneer Venus configuration. 

1.2 Formulation of Problem 

An electric field propagating in free space can be written as 

UQ(r) = Ao(r) e (1) 

where Aq is the amplitude and <f>o the phase. If the medium is turbulent, 
the electric field becomes 
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u(r, t) 


( 2 ) 


A(r,t) e 


t) 


Using Rytov’s method to solve Maxwell’s equations, it can be shown (Ref. 4) 
that 

u(r,t) = Uo(r) 


(3) 


where - x + iS , and X = £n is the log-amplitude fluctuations while 
S = 4» - is the phase fluctuation. Ishimaru (Ref. 3) has shown that if 
B^g is the cross correlation between the log-amplitude and phase fluctuations, 
then the temporal frequency spectrum W^g of the cross correlation between the 
log-amplitude and phase fluctuations defined by 




/ 


^ -imT , 

B _ e dT 
XS 


(4) 


where f is the fluctuation frequency and oj = 27rf, is given by 




/ Kd. 1C 
./V t(Kv)^ - 




(5) 


where k is the spatial wave number, spectrum for the index of 

refraction fluctuations, v is the wind velocity transverse to the line-of- 
sight path, k is the free space wave number and L’ is the propagation distance 
For spherical wave propagation 
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where 


h 


exp 


- 1 




* denotes complex conjugate, n is the coordinate along the line-of-sight path 
originating at the transmitter and Im[ ] means imaginary part of [ J . 

Another quantity we wish to calculate is the coherence coh of log- 

Ao 

amplitude and phase fluctuations (Ref. 5). 


coh^g = p^(£) 


ks^^)r 

W^(f) WgCf) 


where W^(f) is the frequency spectrum for the log-amplitude fluctuations and 
Wg(f) is the frequency spectrum for the phase fluctuations. In many analyses 
the square root of the coherence p(f) is used since it is more analogous to 
the familiar coefficient of correlation (Ref. 6). When p (f ) = 1 the log- 
amplitude and phase fluctuations are completely coherent at f , while if 
p (f ) = 0, they are Independent at f. 


1.3 Application to Pioneer Venus Entry Mission 

The Pioneer Venus configuration is depicted in Figure 1. As discussed 
in our previous report (Ref. 1), we will assume that turbulence in the 
Venusian atmosphere is homogeneous for altitudes lower than 55 km and non- 
existent for higher altitudes. From our study based on Mariner 5 (Ref. 7), 
we have found that the intensity of turbulence actually varies with altitude. 
The effect of this variation on the gross features of the frequency spectra 
should be small, so that the results obtained in this report are representative 
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Another point worth mentioning is that for the Pioneer Venus entry mission, 
the Fresnel-zone size (A is the wavelength and L is defined in Figure 1) 

is smaller than the outer scale size of turbulence Lq (Ref. 1) . Thus, the 
effects of finite extent and inhomogeneity of turbulence in the direction 
transverse to the line-of-sight path that need be considered for the cross 
spectra in a fly-by mission (Ref. 8), can be neglected in the case of Pioneer 
Venus. As in our previous study (Ref. 1), we will assume that 'I'jj(k) is 
represented by Kolmogorov’s spectrum of the form proposed by von Karman, i.e.. 


/ 1 

~ 0.033 c^ /k^ + ~2 y 


( 8 ) 


where c^, the structure constant for the refractive index fluctuations, 
characterizes the intensity of turbulence and Lq is the outer scale size 
of turbulence. 

Let us now apply (5) and (6) to the Pioneer Venus configuration shown 
in Figure 1. We note that L” = L’ + L and since L” >> L, n(L-n)/L 
In (6) becomes simply n. We then integrate (6) and obtain 


Also, (5) can be rewritten 




v 


00 




d<’ 


( 10 ) 


In the next section we will evaluate (10) and subsequently (7) . 
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1.4 Results and Discussion 


Using the second solution of the confluent differential equation 
(Ref. 8), (10) along with (9) can be expressed in the following closed 
form: 


(f) 


0.033 4 ^ rg) 


iil 

\ 6 / 


Re j^e^“ ~ ^ I “ 


where a = -p- > 


( 11 ) 


V /k~ 

= > Rs[ ] means real part of [ ], F is the Gamma 

function and ij^(a, c; z) is the confluent hyper geometric function which is 
independent of the Kummer function. Let us examine the asymptotic forms of 
(10) when f 0 and f ->■ <». 

For the case f ^0, we write ip(a, c; z) in terms of the Kummer function 
i}> (a, c; z) (Ref. 9) 


i()(a, c; z) 


r(i-c) 

r(a-c+i) 


+ 4.(a-c + l, 2-c; z) 


( 12 ) 


and then use the small argument approximation for <|) (a, c; z) (Ref. 9) 


(j)(a, c; z) 


1 + - 2 
c 


(13) 
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We finally obtain 


f 0 


0.033 4 


47t2 7/3 ,2/3 

L k 

V 



32 


-2/3 

a 


(14) 


For convenience of plotting we will normalize (14) to W^(0) as given in 
Ref. 3, i. e. , 

W^(0) = 

We then have 

”xs<°> 

Wx(0) 

For the case f ->“00, we employ the large argument approximation for 
i|)(a, c; z) in (11) (Ref, 9): 

ip(a, c; z) = z ^ (17) 

and again normalizing to W^(0) we obtain 

XS^ ^ 

W^(0) 
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So far we only have the asymptotic solutions normalized to W^(0). 

To calculate the intermediate values we used the Univac 1108 digital computer 
to integrate (5). The method employed was the trapezoidal rule together with 
Romberg's extrapolation method (Ref. 10). Shown in Figure 2 are the results 
of W^g(f) normalized to W^(0) and plotted as a function of a for various 
values of 3 = * These values of 3 correspond to = 50, 100 and ® m 

for L = 55 km at f = 2300 MHz. As can be seen the fluctuations in all cases 
occur primarily at frequencies less than 2fj_, For L = 55 km, and assuming 
v = 50 m/sec, 2fj^ = 0.47 Hz. 

Calculations for the square root of coherence p have also been made and 
the results are plotted in Figure 3. The coherence results show that there is 
no correlation between the log-amplitude and phase fluctuations for f > 

When f < 2 f_L there is some correlation with the peak of 'vO.65 occurring 
at f 

The results discussed so far assume that the sub-earth direction 
corresponds to the zenith of Venus (Figure 1) . If the sub-earth direction 
makes an angle 0 with the zenith of Venus, then L in (10) must be replaced 
by L sec 0, while v by v cos 9. 

1.5 Conclusions 

We have studied the cross correlation between the log-amplitude and 
phase fluctuations of the Pioneer Venus entry probe radio link. It is seen 
that for fluctuation frequencies above 1 Hz or so there is no correlation. 

For frequencies below, the log-amplitude and phase fluctuations are weakly 
correlated with the peak in p occurring at f f and having a value close 
to 0.65. 
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II. The Study of Small-Scale Turbulence in the Atmosphere 
of Venus Using the Mariner 5 Phase Fluctuations 


2.1 Introduction 

As discussed in our previous report (Ref. 1), prediction of the effects 
of turbulence in the Venus atmosphere on the radio link of Pioneer Venus entry 
probes relies heavily on radio measurements of the turbulence made by previous 
space probes. We have studied the amplitude fluctuations of the S-band signal 
received from Mariner 5 when it flew behind Venus. The results yielded the most 
comprehensive measurements of the small-scale turbulence characteristics of the 
Venus atmosphere for altitudes higher than 35 km (Ref, 7). In this report, we 
wish to examine the possibility of deriving similar information based on the 
Mariner 5 phase fluctuations. In Sections 2.2 to 2.5 we discuss the processing 
steps for extracting the turbulence- induced phase fluctuations from the Mariner 5 
radio data, while the results are discussed in Section 2.6. 
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2«2 Extraction of Phase Data from the Received Signal 


The signal transmitted from the Mariner 5 spacecraft through the 
atmosphere of Venus is received on earth at a later time due to delays 
which can be attributed to (1) the propagation time required for the signal 
to propagate from the spacecraft to the earth, (2) the delay caused by the 
"average”, or non- turbulent, Venusian atmosphere, and (3) deviations from 
this average caused by turbulence and other atmospheric non— homogeneities* 

It is this last component which is to be analyzed. 

The delay which is to be measured is manifested in the received signal 
as a phase modulation of the received carrier. The receiving and processing 
steps to which the signal is subjected are designed to extract this phase 
modulation, and then to identify which of its components are due to turbulence. 
Then these components are extracted and spectrum analyzed. 


2.3 Phase-Locked-Loop Processing 

The signal v(t) from the receiving and pre-processing operations 
described in Appendix A of Reference 1 is 


v(t) = A(t) vT sin {2Tr[(A/2)t - fsc<in(t)] + + n(t) (19) 


where 

A(t) 

A 

f sc 


rras received signal amplitude (time varying), 
final signal bandwidth processed (Hz), 
carrier frequency transmitted from spacecraft. 
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d^(t) = delay caused by non-homogeneities in the Venusian 

atmosphere and other uncertainties (time varying), 

= arbitrary constant phase offset (radians), 
n(t) = receiver noise added to the received signal. 

In order to quantitatively determine the effects of planetary turbulence, it is 
necessary to extract (or estimate) the amplitude A(t) [as described in Ref. 1] 
and the phase fsc^n(t) from the signal v(t). This is accomplished by computer 
processing which extracts these parameters with a digital phase-locked-loop 
(PLL) . 

It is important to note that the delay component d^j(t) consists not only 
of the delay caused by non-homogeneities in the Venusian atmosphere, but also 
other uncertainties. These other uncertainties include (1) error in the assumed 
transmitted frequency fgc ho oscillator drift, (2) error in the assumed 
trajectory delay component, and (3) error in the assumed average atmosphere 
delay component (the so-called ’’spline-curve fit error"). These errors must 
be recognized and removed in subsequent processing, and the only basis presently 
available for doing this is the assumption that these errors are more slowly 
varying functions of time than is the turbulence component, so that they can be 
filtered out, i.e. , separation will be done in the spectral domain. 


2.3.1 Phase Extraction 

Here we will not give a complete description of the workings of the 
phase-locked-loop, since several good references on this subject are readily 
available. Instead, we will simply note that the PLL is a non-linear phase- 
tracking device which can normally (when it is tracking) be well approximated 
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as a device whose input is of the form 


x(t) = k\/2 sin 1 2 7T f^j,Q t + 6l(t)| + n(t) 


( 20 ) 


where fyco quiescent frequency of the PLL's VCO, and which extracts the 

phase 01 (t) from the signal. In the process of doing this, the PLL passes the 
phase 01 and the noise n(t) (scaled by division by the amplitude A) through a 
lowpass filter to produce 02, a noisy estimate of the phase 0i, In other words, 
the output phase estimate 02 (t) is that which would be obtained by passing the 
input phase deviation 8i(t) through a lowpass filter, plus an error which can 
be thought of as coming from passing white noise through the same lowpass filter. 

When using the PLL to estimate the signal delay component due to 
atmospheric non-homogeneities, the VCO quiescent frequency is set to the 
center of the band at A/2 so that 


01 (t) = - 2 IT fscdn + kl • 


( 21 ) 


The phase estimate 02 (t) is then simply noise plus a term proportional to the 
desired delay component d^ lowpass filtered. 

The transfer function of the equivalent lowpass filter of the linearized 
second-order PLL is (Reference 11, Eq. 5-19) 



( 22 ) 


1 + (2C/Wn)s 

1 + (2c/Wn)s + s2/wj 
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where 


Bl = one-sided loop noise bandwidth (Hz), 

R = Tausworthe’s damping parameter, 

Wjj = loop natural radian frequency, 

^ = loop damping factor. 

The magnitude and the squared magnitude of the filter transfer function L(j2Trf) 
are shown in Figures 4 and 5, respectively. 


2.3.2 Phase Errors 


The phase estimate 02(t) produced by the PLL will differ from the desired 
phase ©2 “ ^''T^sc'^n ^1> factors, viz., (1) higher-frequency 

variations In S^Ct) will be filtered out by the lowpass filter L(s) (unless the 
loop bandwidth is made sufficiently large), and (2) the low-frequency compo- 
nents of the noise n(t)/A, which in this case has an essentially flat spectrum, 
will pass through the filter L(s). Thus, it is clear that for the best possible 
phase-estimate accuracy, the loop bandwidth used should be as small as possible, 
consistent with the requirement that it pass all significant variations in the 
input phase 0]^. 

If we let 


(f> = Phase error (rad) 

= 0i(t) - 82(t) 

and 


(23) 


S^(f) = one-sided spectral density of <j> due to noise. 


(24) 
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then for spectrally flat noise, which is typical of the receiver noise present 
at the input to the PLL, with one-sided power spectral density Nq watts/Hz, 
we get 




»0 , 

— |L(j2itf) 

A ^ 


( 25 ) 


o 

where Nq/A is the spectral density of the noise at the input to the equivalent 
lowpass filter L(s), and A is the power in the signal into the PLL (see Ref. 12, 
Sec. 2.8), Thus, the spectral density of the noise-induced phase error will have 
the shape shown in Figure 5, and a level which depends on the signal- to-noise 
ratio (SNR) Nq/A^. 

The mean squared value o| of the phase error is just 


f S^(f) df 

0 



(26) 


In order for the linearized model of the PLL which we have been assuming to be 

2 

valid, this mean squared phase error must be less than about 0.1 (rad) , Thus, 
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for the PLL processing described here the loop bandwidth was chosen small 
enough to satisfy 


2 




0.1 


( 27 ) 


B 


L 


- 10 N 


0 


SNR in 1-Hz Bandwidth 
10 


(28) 


2.4 Spectrum Analysis of the Phase Estimates 

Once the phase as a function of time has been estimated by PLL tracking, 
it is desired to separate this estimate into three components; (1) the fluc- 
tuations due to the Venusian atmospheric turbulence, (2) the (slowly varying) 
errors due to spacecraft oscillator drift, trajectory errors, and errors in 
the assumed average atmosphere, and (3) the variations due to noise at the 
receiver input. It was expected that the separation could be made in the 
spectral domain by the following two steps. First, the phase estimates are 
"detrended" by passing them through a high-pass digital filter which removes 
all dc and slowly- varying components from them. Then, the high-pass-filtered 
signal is spectrum analyzed using an FFT algorithm, and theoretical results 
are used to identify the noise components. 

2.4.1 High-Pass Filtering 

Since spectral analysis is a special case of bandpass filtering, it is 
not always necessary to precede the spectral analysis step by high— pass 
filtering or "detrending". However, when this is not done and a large dc 
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or "trend" component is present, the sidelobe response of the bandpass 
filters of the spectral analysis will respond to the very large dc component, 
and the high-frequency portion of the spectrum will be distorted. For the 
phase data spectrum analyzed here, a Butterworth highpass filter of order 4 
to 8 was used when the high-frequency portion of the spectrum was being 
studied, and no filter was used when the low-frequency portion of the 
spectrum was being studied. This was so that when trying to separate the 
residual delay errors from the low-frequency turbulence components in the 
delay d^, there would not be any filter artifacts near its cutoff frequency 
to degrade the spectra. 

2«4«2 Spectrum Analysis 

The samples of phase, with or without highpass filtering, were processed 
through a fast Fourier transform (FFT) routine to estimate the spectrum of the 
signal’s phase. This was a conventional FFT routine augmented by two special 
features: a Hamming window and zero fill. The purpose of the Hamming window 

is to reduce the spurious "sidelobe" response of strong spectral components 
(with a corresponding loss in frequency resolution), and the purpose of using 
zero fill is simply to evaluate the Fourier transform of the phase estimates 
at closer-spaced frequencies than would otherwise be obtained. 

2.4.3 Spectrum Normalization 

For direct comparison wifh the theoretical results it is necessary to 
maintain the correct spectral level at the output of the spectrum-analysis 
processing steps. This means that the effect of the Hamming window must be 
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accounted for, as well as that care must be used in determining the appropriate 
normalization at the other steps. 

The Hamming window is accounted for by computing the mean-squared value 
of the phase data both before and after the Hamming window is applied. Then 
the spectral levels produced by the FFT routine are increased by the ratio of 
these mean-squared values. This assures that the area under the spectral curve 
equals the energy in the data, as is required by Parseval^s theorem. 

2.4.4 Frequency Averaging 

The final step in producing the experimental spectra is averaging over 
frequency. This step is required due to the fact that we are estimating the 
power spectrum of a random process (the turbulence- induced phase fluctuations) 
from one sample function of that random process, and estimates obtained this 
way are very noisy. For example, the standard deviations of such estimates 
are approximately equal to the mean spectral level. By averaging several 
adjacent frequency points together, the variances of the estimates are reduced, 
at the expense of reduced frequency resolution. 

The averaging done here is a "running-point” average, i.e., the spectrum 
at every frequency point is replaced by 

(m— 1) / 2 

^ ^ G(f + tAf) (29) 

where 

m = number of points (assumed odd) in running-point average, 

Af = frequency spacing from FFT, 

G(f) = spectrum from FFT at frequency f. 
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2,5 Interpretation of Results 


The last step in the processing consists of using the theoretical 

results of this report to separate the turbulence-induced fluctuations from 

the noise-induced fluctuations. The noise is recognized [see Equation (25)] 

2 

by having the spectral shape shown in Figure 5 and a level given by Nq/A . 
The part of the experimental spectrum which deviates from this noise part is 
then interpreted in terms of turbulence- induced phenomena. 


2,6 Results and Discussion 

Shown in Figure 6 is the time history of the Mariner 5 signal level during 
entrance occultation. The altitudes corresponding to the lowest penetration of 
the propagation path are indicated at the top of the figure. As discussed in 
Ref, 7, turbulence was found in regions A and C, which correspond to the vicin- 
ities of 60 and 45 km, respectively. Region E is of interest also because this 
pre-occultation period indicates the "system” noise. Shown in Figure 7 are the 
frequency spectra of the log-amplitude fluctuations which were obtained in our 
previous study. It is clear that the spectral density levels for regions A and 
C are substantially above the spectral density levels of region E and the noise 
levels corresponding to regions A and C as well. The frequency spectra obtained 
for the phase fluctuations in regions A, C and E are depicted in Figure 8. The 
frequency response of the eighth order Butterworth high-pass filter is shown in 
Figure 9. The cut-off frequency is 0.1 Hz. As discussed in Section 2,5, the 
corresponding noise levels are given by Nq/A^. These are not shown in Figure 8 
because they are lower than the range of spectral density levels plotted. The 
theoretical plot is the theoretical frequency spectrum for phase fluctuations 
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for region A. This has been calculated using the analytical results given 
in Ref. 13 and the turbulence characteristics derived from the log-amplitude 
fluctuations (Ref. 7). Contrary to the log-amplitude fluctuations, the 
frequency spectra for the phase fluctuations for regions A, C and E show no 
substantial differences. Moreover, the spectral density levels are signifi- 
cantly higher than the theoretically predicted ones. These results indicate 
that the errors due to oscillator drift, assumed trajectory delay component 
and spline-curve fit exceed the turbulence induced fluctuations and have not 
been completely removed by our presently developed processing techniques. 
Looking at it from a total power point of view, the standard deviation of 
phase fluctuations due to the errors mentioned above exceeds the standard 
deviation of phase fluctuations due to turbulence which is estimated to be 
approximately 6 degrees. 

It should be mentioned that a study of Venus atmospheric turbulence 
based on the Mariner 10 phase fluctuations is currently being carried out. 
Because of the availability of a coherent S/X link on Mariner 10, the 
uncertainties present in the Mariner 5 data will hopefully be removed by 
looking at the phase differences. 


2.7 Conclusions 

We have studied the Mariner 5 phase fluctuations with the objective 
of inferring the turbulence characteristics of the Venus atmosphere. It 
is seen that with the data processing techniques developed and currently 
available, the phase error due to oscillator drift, assumed trajectory 
delay and spline-curve fit exceed the turbulence induced fluctuations in 
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the phase. It is therefore not possible to infer the turbulence characteristics 
from the Mariner 5 phase fluctuations. Our study of Mariner 10 should, however, 
shed more light on this problem, since these errors can be eliminated by looking 
at the S/X phase difference. 
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Fig. 2, Normalized cross power 
spectra of log-amplitude and phase 
fluctuations 


Fig. 3. Square root of coherence 
of log -amplitude and phase 
fluctuations 
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Fig. 6. Time history of the signal level for Mariner 5 entrance 
occultation (from Ref. 7) 
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Fig. 8. Frequency spectra of the 
Mariner 5 phase fluctuations for 
regions A, C and E of entrance 
occulta tion 



Fig» 9. Frequency response of 
eighth-order Butterworth high pass 
filter. 
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